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m ' 1. INTRODUCTION 

Our main goals have been to calculate photonic bands on various periodic dielectric 
lattices, to clarify some controversy regarding existence of a photonic band gap for fee 
lattice of dielectric spheres and finally treat impurities in a photonic crystal 0. 
Motivated by the search for a photonic band gap || we have tried to adapt one of the 
standard methods of electron band theory - namely the scalar bulk Kohn-Korringa- 
Rostocker (KKR) method [§, || to the form appropriate for photons @, ||, . We recall 
that in the case of photons the role of a (periodic) potential plays v(r) = e(r) — e G . 
Here e(r) is the dielectricity of a medium and e Q is its host value which is assumed 
to be uniform and homogeneous || |(|. A dielectric "atom" V s is called a connected 
region where v(r) ^ 0. The reason of our choice was that the KKR method proceeds 
analytically as far as possible and enables to go beyond nearly-free photon and plane 
wave approximations. The latter were cast into doubts for discontinuous potentials [jl], 
and by a rigorous proof on the existence of finite number of gaps in two and higher 
dimensions for the periodic Schrodinger operator ||. The KKR starts with an integral 
equation which is after expansion in a suitable basis transformed into an algebraic one. 
A band structure then follows from the conditions of solvability of the algebraic equation 
which is vanishing of a determinant of a matrix (see below) which determines dispersion 
relation and eventually photonic bands. The electronic KKR method is known to lead to 
a very compact scheme if the perturbing periodic potential v(r) is spherically symmetric 
within inscribed spheres and zero (constant) elsewhere ||, ||. In the case of electrons 
already p-wave approximation gives agreement within 2% with experiment. 

Another distinguished feature of the method is a separation of pure geometrical and 
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scattering properties of a medium. Geometrical properties are encoded in geometrical 
structure constants Ai m -if m f, Ai> m / ; i m = Ai m .i> m /*, characteristic for the lattice under 
consideration. They are functions of energy a (for photons a = uj ^/e^, where uj is 
frequency) and the Bloch momentum k. Once the structure constants are known then 
all that is needed in the case of lattice of identical scatterers is to know scattering 
properties (phase shifts) of a single scatterer. 

We shall try to elucidate you a derivation of the photonic KKR along the lines of 
Q and from the multiple-scattering theory (MST) ||. The latter is a formal method 
for calculation of the spectrum in a (photonic) crystal under the presence of impurities. 
In fact it gives a formal solution to the problem of finding the spectrum in the case 
of a collection of arbitrary distributed scatterers of different shape and of different 
scattering properties. In order to simplify our discussion dielectric medium is assumed 
to be magnetically isotropic and from now on we shall set /i(r) = 1 as well as c = 1. On 
discontinuities S of e(r) side limits of e(r) and its derivatives as well as limiting values 
of fields and their derivatives are assumed to be well defined, too. 



2. INTEGRAL EQUATIONS 

In course of their derivation we have found that there still persists some confusion. 
This originates from the fact that unlike the scalar case in the case of photons fields can 
change dis continuously across discontinuities of electric and magnetic permeabilities and 
fields are different on different sides of discontinuities. This is the essential difference 
between the Schrodinger and the Maxwell equations^. Therefore, in the latter case one 
has to carefully distinguish between "inward' and " outward 1 formulations, i.e., whether 
in a given surface integral inward or outward limits of fields and their derivatives at the 
atom boundary are taken [fj. |(| [7J . The inward formalism is basically that of KKR ||] 
while outward formalism is originally due to Morse ||. In particular we have found that 



the photonic KKR method proposed in [T(J does not serve the purpose. There fields in 



integral equations are outward limits of fields with respect to the atom boundary while 
a scalar product with inward limits of fields is taken. 

The outward integral equation for E(r) can be written as follows |J, 

E(r) = / [E(r') (dS' ■ V)G a (r, r') - G a (r, r') (dS' • V')E(r')]. (1) 

•>dV s+ 

Here G a (r,r') is the Green function of the Helmholtz equation in a region Q with 
suitable boundary conditions imposed, defined by 

(V 2 + a 2 )a(r, t') = 5(t-t'). (2) 

The same equation also holds for B(r). Thus, formally both B(r) and E(r) satisfy the 
same integral equation ([I]) outside a dielectric atom V s . What makes the difference 
between them are matching conditions across a boundary of the atom. 
Corresponding inward integral equation for E(r) is |J 

E(r) = XVs {v) E(r) + / [G CT (r, r') (dS' • V')E(r') - E(r')(dS' • V') G a {r, r')] 

- - I [V'G CT (r, r')Mr') (E(r') • dS') - / G a (r, r') [V • E(r')] dS', (3) 



3 Another principle difference is that the Maxwell equations do not seperate even if e(r) is a sum of 
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Xv a ( r ) being the characteristic function of V s || . 

The equivalence of the two expressions can be established either by taking different 
limits in integral equations for the vector potential A(r) || or directly by looking at 
the matching rules of fields and their derivatives on different sides of E. It can be found 
in almost any textbook that 

(K~E~)\ = e -^±E- = V -^E~, Et = E-. (4) 

However, none of the standard textbooks (as far as we know) tells you how derivatives of 
fields behave on a discontinuity of permeabilities. Our results are that normal derivative 
<9 n E 4 (r) of tangential component of E(r) changes discontinously across E, 

(d n E+ - a„E-)| =(---) V t D n = -V t E~, (5) 

V t D n (r) is continuous and d n Ef = ^ t E~. Rather surprisingly, albeit normal compo- 
nent E n (r) of E(r) is discontinuous, 

(d n E+ - d n E~)\ x = - (j^(d n e) + - ^(d n e)^j , (6) 

and provided {d n e) + = (9 n e)_ = the derivative changes continously |J. In this manner 
singularities appearing due to discontinuities of permeabilities are safely treated. These 
relations explain why the Gauss theorem which requires at least continuous functions 
can be applied to electromagnetic fields despite that fields change discontinuosly. 



3. PHOTONIC KKR METHOD 

In the case of photons the KKR method has been only used within the scalar ap- 



proximation to the Maxwell equations |TT[ . Following the lines of [|J] we have formulated 
a general variational principle for the Maxwell equations which holds for an arbitrary 
shape of the basic "atom" of a dielectric lattice. The photonic analogue of the scalar 
KKR functional is defined to be 

A = u 2 I dhv(r)E*(r) • (e(i-) + u 2 f G (f (r,r>(r')E(r , )dY+ 

+ [ Va^rOfV'-E^rfV). (7) 

The integral over r is well defined and does exist as a bothside limit. Provided e(r) is 
real (no absorbtion), and using well-known hermitian properties of Green's functions, 
d T G a (r, r') = — d r >G a (r, r') and G* (r, r') = G CT (r', r) one can check that variations of A 
with respect to E(r) or E*(r) reproduce correctly the equation for electric field E(r) or 
its complex conjugate E*(r) within V s , respectively. Rather surprisingly, unless e(r) is 
piecewise constant, i.e., the classical muffin-tin potential, there is impossible (at least 
for our variational principle) to write (0) in terms of surface integrals only. Thus, when 
working with the variational KKR functional it is necesary to confine ourselves to the 
case when e(r) is nothing but a real muffin-tin potential. Then in the case of spherically 
symmetric dielectric atoms general formula (|7|) can be simplified further as follows, 

A := lim / dS [d r E*{r) - E*{r)d r ] -U dS' [^E(r') - E(r') d r ,\ GJr, r) 

-- I [V'G a {v,v')\v{v'){E{v') ■ dS')\ . (8) 

e JdVs- 



The final expression (§) resembles the scalar case Q, the only difference being the last 
term. The formula (|7]) is called off-shell while the later (|3|) is on-shell. The reason 
behind this is that the latter formula is derived from (|7|) by using field equations which 
for a (relativistic) free particle means that it is on its mass-shell. Therefore in looking 
for an extremum of (0) one is allowed to use " arbitrary" test functions while in the 
case of (|]) one has to ensure that the test functions are local solutions to the Maxwell 
equations. In the latter case the variational principle selects among the local solutions 
that which satisfies prescribed (Bloch) boundary conditions. 

To find the local solutions is rather difficult. Usually only inward solutions can be 
found. From now on the further treatment along lines in [|j] is straightforward, albeit 
more involved @. To proceed further analytically E(r) and G a (r, r') are expanded into 
spherical harmonics [|T^, [13] after which the integral equations are transformed into 
matrix equations for expansion coefficients. If it is done in (|]) one obtains (in our 
notation) direct KKR method while starting from (|8|) variational KKR method 0, [|. 
The "direct" photonic KKR method has a wider region of applications : complex and 
nonconstant e(r) within atoms, while the variational KKR method has a rather limited 
range of applications. However, being variational it is expected to converge more rapidly 
within its range of application. 



4. MULTIPLE-SCATTERING THEORY 

Standard assumption of MST is that a given medium can be divided into non- 
overlapping spheres V n ' s each of which contains one and only one scatterer. MST 
implicitely assumes that the t-matrix (transition factor ||) or phase shifts for a single 
scatterer are known. Once they are known it is convenient to use the outward integral 
equation which is simpler than the inward one. However, one must not forget that the 
problem of finding inward local solutions still remains - they are necessary to determine 
just the t-matrix! Following a standard derivation || one starts with the on — shell 
Lippmann-Schwinger equation for E(r) 0, 

E (r) = W dS'[d r ,G a (r, r') - G a (r, r') d r ,]E(r'), (9) 

which holds provided r stays inside spheres. E D (r) is an incident wave and n labels 
scatterers (or spheres containing single scatterers, as you want). After expanding fields 
and the Green function into spherical harmonics one arrives at the basic photonic MST 
equations @, 



jA'L' 



0~ijO~LL'0~AA' + i 2^ G ALA » L „t A „ L „ A , L , 
A" L" 



C jA/L >. (10) 



Here indices i, j label different scatterers, G n = 0, and C° AL or Cjal are expansion 
constants of E G (r) or E(r) into spherical harmonics around j-th sphere, respectively. 
In contrast to the scalar case the expansion is more involved due to vectorial nature 
of fields. This is reason for additional index A which labels magnetic (A = M) and 
electric (A = E) multipoles As a result the order of matrix ( |TUD in a given (s-, 
p-,. . .-wave approximation) is doubled with regard to the scalar case. G %3 AL A „ L „ is the 
vector (propagator) Green function. Physically it describes what amount of a particular 
multipole field scattered from the z-th site contributes to the particular multipole field 
incident on the j-th site. 



To write the vector Green function G l AL A , L , explicitely we have defined quantities 
C a in terms of 3j-symbols C a {l'm!lm) [13], 



C a (l - l,m + a) = Vl + lC a (l-l,m + aJ,m), 

C a (l + l,m + a) = -VlC a (l + l,m + a,l,m), (11) 



Tl m + a — r— ^Im+ct) (12) 

'1(1 + 1) 



and 



where T£ are defined by action of spherical components L a of orbital angular momentum 
on spherical harmonics Y L , L a Y L = T^ +ct Yi m+a [||, [7J. Radial functions which are 
proportional to ji(o~ri)' s (the spherical Bessel functions) are normalized to unity on the 
boundary of the z-th sphere. This simply amounts to redefinition of C° ML and Ciml- 
We also redefine G l EL EL , by dividing factor a 2 to which Jel s are normalized (they have 
different dimension than 3ml' s because of rotation operation applied). Afterwords the 
resulting expresions can be written in the following compact form 

1 

GmL,ML' = i E 9l'm> '+a;lm+a^hn+a^Vm' ' +ai (13) 
a=-l 

1 1 

G ML,EL> = ° E E 9 l l? +p >,m'+a,lm+<xC a ( 1 ' + P' > m ' + a )^m+a> (14) 
p'=-l a=—l 

1 1 1 



GeL,ML' ~ ~ E E 9l',m'+a,l+p,m+aC a (l + P, m + Ot)T^ m , +a1 (15) 
°" p=-la=-l 

G%L,EL> = « E E E^,m^l4*n^^(i+P^ + «)C m (i / +^ "»' + «)• (16) 



a=— 1 p=— i »'= 



Note that is not Hermitian. This is a consequence of using the f-matrix which 

also is not Hermitian j|, [|. For spherically symmetric scatterers the ^-matrix is diagonal 



and for homogeneous spheres it is explicitely known as the Mie solution [|T4 |. 

Provided scatterers are identical and arranged in a periodic way one can take 
Fourier transform with regard to the Bloch momentum k. The condition of existence 
of a solution to (|10D for C° AL = 0, 



det 



5lL"5aA' + i E GAL,A"L"(k)tA"L",A'L' 



0, (17) 



then gives the photonic KKR equation. 

We have noted an attempt to derive MST for photons HT5|| . However one can find 



obvious difference : our expression is much more symmetric while in [To] there is no 
summation over p, p' which is necessary here since electric multipoles of order / have 
nonzero matrix elements with Yjy only for V = I ± 1. Moreover, in our case vector 
structure constant are determined via 3j-symbols and T" coefficients in contrast to the 
Gaunt numbers in |T5| . 

Presented results use as much as possible existing electron structure constants 
and hence allow for a direct numerical application. They allow to calculate photonic 
bands and to treat impurities as well. Despite that a photonic band calculation is more 
involved due to its vectorial nature there is one small advantage with regard to the 



Schrodinger equation : the Maxwell equations are conformally invariant and therefore 
a scale of a dielectric lattice can be chosen arbitrary. By taking the scale identical to 
that of electrons the scalar structure constants can be used without any changes. For 
more complete treatment of the above problems together with numerical results see 0. 

Last but not least some new phenomenon in a photonic crystal which has not been 
discussed yet and which may be of some importance : particles propagating through 
such medium can emit Tcherenkov radiation at rather small velocity. 
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